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: I. INTRODUCTION 

p i' Exact solutions played very important roles in the development of physics. The exact 

solutions of the Schrodinger equation for a hydrogen atom and for a harmonic oscillator 
^ . laid a foundation of quantum mechanics. The next simplest atom is the helium atom, 
^ ! which is a typical three-body problem and has not been well-solved both in classical 
lO I mechanics [1] and in quantum mechanics [2]. The Faddeev equations [3] are popularly 
^ ■ used in both scattering processes [4] and bound state calculations [5] for the quantum 
O ■ three-body problem. However, only a few analytically solvable examples were found [6]. 
^ . For the helium atom, variational methods can achieve 9- to 12-place precision for energy 
! values with a few hundred or thousand variational parameters [7-10]. There are many 
O I reasons to prefer the direct solutions of the three-body Schrodinger equation over the 
c/3 ■ variational one, for example, the analytic structure of the variational wavefunction is 
chosen arbitrarily. The accurate direct solution of the three-body Schrodinger equation 
P^. with the separated center-of-mass motion has been sought based on different numerical 
^ ! methods, such as the finite difference [11], finite element [12], complex coordinate rotation 

' l'^'^]^ hyperspherical coordinate [14], and hyperspherical harmonic [15-17] methods. 

^ ' In the hyperspherical harmonic method [15-17], the six Jacobi coordinates, after sep- 

aration of the center-of-mass motion, are separated to one hyperradial variable p and five 
hyperangular variables fl. The wave function is presented as a sum of products of hyper- 
radial function and the hyperspherical harmonic function, depending on Q. Since three 
degrees of freedom for the rotation of the system are not separated completely from the 
internal ones, there is the huge degeneracy of the hyperspherical basis. The interactions 
in the three-body problem is not hyperspherically symmetric so that the matrix elements 
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of the potential have to be calculated between different hyperspherical harmonic states. 
Another difficulty in the practical calculations by the hyperspherical harmonic method 
is the slow convergence of the series. The convergence was fastened by decomposing the 
wavefunction ^ = x0, where x is chosen to take into account the singularities of the po- 
tential and the clustering properties of the wavefunction, and is the part to be expanded 
in hypersphericals [15]. 

In an unpublished paper [18], Hsiang and Hsiang devised a principle method to reduce 
the Schrodinger equation of three-body problem to a system of linear algebraic equations. 
Although their method was incomplete, it contained a good idea for simplifying the quan- 
tum three-body problem. Developing their idea, we suggest an improved hyperspherical 
harmonic method to separate three rotational degrees of freedom completely from the 
internal ones, and the decomposition of the function for fastening the convergence of the 
series in the hyperspherical harmonic method [15-17] is still effective in the improved 
method. 

First, after the separation of the center-of-mass motion and the whole rotation, the 
three-body Schrodinger equation with any given angular momentum and parity is reduced 
to the coupled partial differential equations, depending upon only three variables, which 
are invariant in the rotation of the system. The coupled equations are closed because 
the potential depends only upon those three variables due to its symmetry in the SO (3) 
rotation and the space inversion. Then, we expand the solutions of the coupled equations 
by the orthogonal bases of functions to reduce the coupled partial differential equations to 
the coupled ordinary differential equations. The bases of functions are similar to those in 
the S wave case of the hyperspherical harmonic method [16], but in our method they are 
effective for any given angular momentum. The coupled ordinary differential equations 
are similar to the hyperradial equations for the hyperradius in the hyperspherical har- 
monic method [16], but there is an important difference between them. The differential 
operators in each ordinary differential equation of our coupled equations are same so that 
the equations can be reduced to a system of linear algebraic equations by the spectral de- 
composition of the hyperradial function. This improved hyperspherical harmonic method 
is effective for the three-body problems of both identical and non-identical particles, and 
effective for any pair potential, which depends only upon the distance of each pair of 
particles. 

The improved hyperspherical harmonic method greatly decreases the calculating time, 
because to solve a system of linear algebraic equations even with the equation number 
more than 1000 is faster than to solve a system of ordinary differential equations with the 
equation number less than ten. Furthermore, all the technique used in the hyperspherical 
harmonic method to fasten the convergence of the series can also be used in the improved 
hyperspherical harmonic method. 

In this paper we will demonstrate the improved hyperspherical harmonic method in 
some detail. In Sec. II we introduce our notations and write the Schrodinger equation 
for a three-body problem in the center-of-mass frame by the Jacobi coordinates. Due 
to the 5*0 (3) symmetry of the Schrodinger equation, in Sec. Ill we will explain how 
to separate the rotational degrees of freedom completely from the internal ones, and 
derive the formulas for reducing the three-body Schrodinger equation with any given 
angular momentum and parity to a system of partial differential equations. In Sec. IV 
we introduce three new variables through the complex vector coordinates [19]. By a 
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complete set of normalized and orthogonal bases of functions of those new variables, the 
Schrodinger equation for any angular momentum and parity is further reduced to a system 
of linear algebraic equations with some coefficients to be calculated. In terms of some 
identities, proved in Appendix, the alternating series contained in the formulas for the 
coefficients are calculated analytically. As an important example in quantum mechanics, 
the problem of a helium atom is studied in some detail. We solve the system of linear 
algebraic equations for the lowest-energy states of the S wave and the P wave in both the 
parahclium and the orthohelium in Sec. V, respectively. Wc also calculate the energies of 
the ground state of some helium-like ions there. The solution of the Schrodinger equation 
for the three-body problem is called exact in the meaning that it is expressed as an infinite 
series in terms of the orthogonal bases of functions, and can be calculated in any precision 
if enough terms in the series are included. Due to our limited condition in computer, in 
the present paper we still use the truncation method of the series and do not fasten the 
convergence of the series by decomposition of the wavefunction, so that the calculating 
precision is in the thousandths for the energy of the ground state of the helium atom 
and in the hundredths for the remaining calculating results. The calculating precision 
will be improved elsewhere. The main purpose of this paper is to present the improved 
hyperspherical harmonic method. Some discussions are given in Sec. VI. 



II. THE SCHRODINGER EQUATION 



Denote by rj and by Mj, j — 1,2,3, the position vectors and the masses of three 
particles in a three-body problem, respectively. The relative masses are mj = Mj/M, 
where M is the total mass, M — J2Mj. The Schrodinger equation for the three-body 
problem is 

-{h^/2M) + = E<if, (1) 
A = Em-^A,^., (2) 
where A^^. is the Laplace operator with respect to the position vector r^- 

A -^ + ^ + ^ 
drji drj2 dvj^^ 

and y is a pair potential, depending only upon the distance of each pair of particles. For 
definiteness, we discuss the Coulomb potential: 

y^Z^^Z^^Z^^ (4) 

where ZjC denotes the electric charge of the jth particle. 

Now, we replace the position vectors r, by the Jacobi coordinates R,-: 



Ri = miri m2r2 + msFa, Ra^x-^J ■ Ri, R3 = y, (5) 

V 1712+1713 



mi 17121713 

x = ri, y=J ra-rg), (6) 

\ 1712 + 1713 V 777,2 + 777,3 
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where Ri describes the position of the center-of-mass and vanishes in the center-of-mass 
frame, x and y describe the position vectors of three particles in the center-of-mass frame 
for two special cases, respectively: 

i) The second particle coincides with the third particle, 



m2 + m3 mi 
ri = -J X, rs = rs = J x, }^ m^r ■ = x . 7 

V mi \^ 1712 + 1713 j^l 

ii) The first particle is located at the origin (the center-of-mass). 
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/ 77227723 , ^ „ „ 

ri = 0, m2r2 = -marg = J -— y, Z^mjr.^y. (8) 

\ 1712 + 1713 p{ ' ■' 

A straightforward calculation by replacement of variables shows that the Laplace op- 
erator in Eq. (2) and the angular momentum operator L are directly expressed with 
respect to R^: 

3 3 

i=i i=i 
In the center-of-mass frame, they further reduce to those with respect to the position 
vectors x and y: 

A = Ay, L = -ihyi x Vx - x Vy (9) 



III. EIGENFUNCTIONS OF THE ANGULAR MOMENTUM 

The Schrodinger equation (1) is spherically symmetric so that its solution can be 
factorized into a product of an eigenfunction of the angular momentum L and a "radial" 
function, which only depends upon three variables, invariant in the rotation of the system: 

6 = X ■ X, ^2 = y • y, 6 = X • y. (10) 

We call them the "radial" variables in this paper. 

For a particle moving in a central field, the eigenfunction of the angular momentum 
is the spherical harmonic function Y^iO^tp). What is the generalization of the spherical 
harmonic function to the three-body problem? 

A natural idea for generalization is to introduce the angular variables, for example, the 
Euler angles. In this way, as discussed by Wigner (see p. 214 in [20]), the eigenfunction of 
the angular momentum is the representation matrix of the 5*0(3) group, D^^^^, («,/?, 7). 
However, the Schrodinger equation with respect to the Euler angles is singular. In the 
hyperspherical harmonic method [16] the spherical harmonic function is generalized to 
the hyperspherical harmonic function, which depends upon five angular variables Vt — 
[a. 9,j.^ 9y^ Lpy), where 6^ {9y) and ip^ (py) arc the angular variables of x (y). Obviously, 
the rotational degrees of freedom are not separated from the internal ones. Is it necessary 
to introduce the angular variables explicitly? 
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In order to avoid the angular variables, let us study the other properties of the spherical 
harmonic function, in addition to the eigenfunction of the angular momentum. As is well 
known, y^i'x.) = r^Yj^{6,Lp), where {r,6,ip) are the spherical coordinates for the position 
vector X, is a homogeneous polynomial of degree i with respect to the components of x, 
and satisfies the Laplace equation as well as the eigen-equation for the angular momentum, 
but does not contain the angular variables and the factor explicitly. The number of 
linearly independent homogeneous polynomials of degree £ with respect to the components 

of xis N{ey. 

e 

N{€) = - s + 1) = + 1)(^ + 2)/2. 

s=0 

Removing those homogeneous polynomials containing the factor = x • x, we obtain that 
the number of linearly independent homogeneous polynomials of degree not containing 
the factor r^, is 

N{e)-N{l-2)^2l + l. (11) 

It is nothing but the number of 3^4 (^) with the same angular momentum I. 

Now, for the three-body problem, there are three radial variables in the center- 
of-mass frame. 3^m(^)3^ni^''(y) ^^e linearly independent homogeneous polynomials of 
degree £ with respect to the components of the position vectors x and y, not containing 
the factors and ^2- The number of those polynomials is M(£): 

M{€) = ^ {2q + 1)(2£ - 2g + 1) = (£ + l)(2f + U + 3)/3. 

Removing those polynomials containing the factor ^3, we obtain that the number of lin- 
early independent homogeneous polynomials of degree £, not containing the factors ^j, 
is 

M{i) - M{e - 2) = 4f + 2. (12) 

According to the theory of angular momentum [20], the polynomials 3^m(^)3^nt'''(y) 
can be combined to be the eigenfunctions of by the Clebsch-Gordan coefficients. The 
number of the combinations with the angular momentum i and {i — 1), which do not 
containing the factors {3, is 

{2£ + 1) (£ + 1) + (2£ - 1) (£ - 1) = 4f + 2. 

It coincides with the number given in Eq. (12). In other words, the eigenfunctions of the 
angular momentum with the eigenvalue £{£+!), not containing the factors ^j, are those 
homogeneous polynomials of degree i or degree {£ + 1). Let us introduce a parameter 
A = or 1 to identify them. Due to the property of the spherical symmetry, what we 
need is to write the eigenfunctions of with the largest eigenvalue of L3, where the 
normalization factors do not matter with us. Denote them by (5g^(x, y) with the degree 
{£ + X), X<q<£. 

A = o: y',{^)yt',{y), 
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Since 3^|(x) ~ {xi + ix^f and J^/_i(x) ~ — \/2l(xi + ix^f ^ X3, we have 
(x, y) = {(g - A)!(£ - + iX2T\y^ + zys)'"^ 

X {(xi+ix2)y3- 2:3(^1 + ^y2)}^ A<?<^, A = 0,1. (14) 

gg'^(x, y) is the common eigenfunction of L^, L3, A^;, Ay, l\xy^ and the parity with the 
eigenvalues ^(£+1), ^, 0, 0, 0, and (—1)^+'^, respectively, where L^, L3 are the total angular 
momentum operators [see Eq. (9)], A^ and Aj, are the Laplace operators respectively 
with respect to the position vectors x and y [see Eq. (3)], and A^j, is defined as 

A _ I I n5) 

dxxdyx dx^dy^ dx-^dy-i 

Their partners with the smaller eigenvalues of L3 can be calculated from them by the 
lowering operator L_ (see Eq. (9) and [20]). 

Now, the solutions to the Schrodinger equation (1), which are the common eigenfunc- 
tions of L^, L3 and the parity with the eigenvalues l{t +!),£, and (—1)^"'"'^, respectively, 
are generally written as 

*a(x, y) = E (x, y), A = 0, 1, (16) 

where ^/^^^ are the functions of the radial variables ^j. Due to conservation of the angular 
momentum and parity, ^^^^(x, y) with different subscripts I and A are separated in the 
Schrodinger equation. Recall that for S wave (^ = 0), go°(x, y) = 1, and ^oo(x, y) = V'o^- 
The wavefunction for S wave only depends upon the radial variables. 

Substituting the wavefunctions (16) into the Schrodinger equation (1) and (9), we 
obtain the radial equations for the radial functions V'^^: 

- ^ AC + 4,^ + 4(. - , + A)^ + 2(, - A)^ + 2(. - ,)-^^-^ 1 



2M i "96 ^ ' 96 ^ ' 96 ^ ' 96 J 

^{E-V)i,'^, \<q<l, A = 0,1. (17) 



IV. RADIAL VARIABLES 

In order to separate the radial variables in the radial equations (17), we define new 
radial variables p, a and /3 through the complex vector coordinates [19]: 

4A;^ (x + iyf = -p'^e~'^ sin a, 

p^2k {6 + 6}'/' , cos a = 2 (66 - 6') (ei + 6)"' , 

tan/? = 26 (6 - Ci)"' , sin/? = 26 [(6 - 6)' + 46']"'^' ' 
k"^ = -2ME/h^, 0<p<oo, 0<q;< 7r/2, - tt < /3 < tt. (18) 



6 



Through replacement of variables, the Laplace operator A in Eq. (9) becomes 

A^——p'— + ^T, 19 
P'^ op op p'^ 

Id (9 1 (9^ 

sin(2a)— + -^ — . (20) 



sin(2Q;) da da sin a d(3 

Denoting by Z„^^(q;,/3) the eigenfunction of T, and expressing it as 

^n,m(«,/3) =e^"^^(sina)l"^lF(C), C = sin^ «, 

we have 

rZ„,^(a,/3) = e^"^^(sina)l"^l 

{^^^^ - 0^ + 4 [(|m| + 1) - (|m| + 2)C] | - \m\{\m\ + 2)| F(C). 
Let F(C) be the hypergeometric function F{—n, n + \m\ + 1, \m\ + 1, (), satisfying: 

|C(1 - 0^ + [{H + 1) - (|m| + 2)C] ^1 F{-n, n + \m\ + 1, |m| + 1, C) 

= —n{n + \m\ + l)F{—n, n + \m\ + 1, \m\ + 1, C)- 

Thus, 

TZn^^{a,p) = -A„^Z„,^(q;,/3), A„„ = (2n + |m|)(2n + |m| +2), 

\ 71 J ^ r!(n — r)!(|m| + rj! 

where the normalization factor is included in front of Zn,m{c(, (^) so that it satisfies 

dp daSmaC0SaZn,mi0i,P)*Zn',m'i0i,P) = 6nn'Smm'- (22) 

-TT JO 

The orthogonal condition is obvious from Eq. (21), and the normalization condition can 
be proved by the identity (A2) in Appendix: 

fTT i-Tr/2 



/TV P'^/ 
dp / da sin a COS a\Zn,m{<^: P)\^ 
-TT Jo ' 



= {2n+\m\ + l)Y^ 



{-iy{n+\ 


\m\ 


+ s)! 


s\{n — s)\{\m\ 


+ s)\{r - 




\m\ 


+ 1) 



X 



{-lY{n+ \m\ 4 


-r)! 




r!(n — r)!( m + 


-r)! 




(— l)"r(r — 


I).. 


• (r — n + 1) 



(r + |m| + l)(r + |m| + 2) • • • (r + |m| + n + 1) 
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where only one term (r = n) in the summation is non-vanishing. As a matter of fact, 
Znm{<^:P) can also be expressed by the Wigner D-f unction [21,16]: 

ZnAo^, P) - (^^^^^) e^'"'^?-r/2)W2) (2«)> when m > 0. (23) 

The Coulomb potential (4) depends only upon the radial variables. It is straightfor- 
ward to show 

|ri-r,r = ^^^i±^^ {l-sinacos(/3-/33)}, (24) 

17111712 ok'^ 

and those formulas obtained by replacing the subscripts (1, 2, 3) cyclically. (3^ is the radial 
variable j3 when the first particle coincides with the second particle, and similar for j3i 
and 132. The calculation results are 

. 2^171117121713 1713- mil7l2 
/3i = TT, Sm/32 = — , COS/32 = ■ 

ms + mim2 ms + mini2 

. -2^/171117121713 1712-71711713 , - 

smPs^ ^ , COS/33 = : (25) 

1712 + lTlim3 m2 + 17li17l3 

By the way, when any two particles coincide with each other, x || y and a = 7r/2. 
Now, taking the Fourier-series expansion for (1 — sin a cos/3)~^/^, we have 

oo 

(l-sinacos/^)"^/^ = ^ T^(sina)e*"'^, 

r=— oo 

TC ^ T (■ \ (4t + 2r)!(sina)^*+'- 

TJsma) = r_r smo; — > — — — t-, rr, rr, 26 

^ ^ ^ ' "to 82*+'^t!(t + r)!(2i + r)!' ^ ' 

Substituting Eqs. (24) and (26) into Eq. (4), we obtain 

-j^—V{p, a, (5)Zn',m'{oi, = Y1 Zn,m{oi, f3)C{\m' - m\)D{n, m, n', m'), (27) 
yoke nm 

D{n, m, n', m') — D{n, —m, n', —m') — D{n', m', n, m) 

— ^2 <i/3 / da sin a cos aZn^rnio^: P)*Tr{sma)e''^^ Zn',m'{<^: P) 

^— ^ J-TT Jo 



[{2n + |m| + l)(2n' + \m'\ + 1)]^^^ ^ 



oo 



g-2t- 


-\^-^'\(At + 2\ 


\m — m'\y. 


t\{t+\ 


\m-m'\)\{2t+\ 


\m — m'\y. 



" {-iy{n+\m\+ry. " (-l)^(n' + |m'| + 5)! 



^ r\{n — ry.{\m\ + ry. ^ s!(n' — s)!(|m'| -|- s)! 

X [t + r + s + l + {\7n\ + \m'\ + \m-m'\)/2]~^, (28) 



C{n) = IJI^ Z2Z3e-^^ + IJ!^ Z3Z,e--^^ + I Z,Z2e-^^. (29) 

y 1712 + ^3 y ITT'S + n7i \j TUi + 7n2 

The formula for D{n, m, n', i7i') contains double alternating series, which can be calculated 
analytically in terms of the identities, given in Appendix. For definiteness, we assume 
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n > n', and first calculate the summation over s by Eq. (A2). Otherwise, we exchange 
{n,m) from {n\m') in D{n,m,n\m'). 



^ {-lY{n' +\m'\+s)\ 



E 



^ s!(n'-s)!(|m'| + s)! 



[i + r + s + l + (|m| + |m'| + |m - m'|)/2] ^ 



(-l)'^'(lm 


+ 6 + r)( m 


+ 6 + r + 1 ) ■ • ■ ( 777 


+ 6 + r + n' - 1) 


( m 


+ a + r)( m 


+ a + r + 1) ■ ■ ■ ( m 


+ a + r + n') 



a = t + 1 + (|m — m'\ - \m\ + |m'|)/2, b — a — n' — |m'|. 

Then, in terms of the identities in Appendix, we are able to express the summation over 
r analytically. 

The remaining differential operators in the radial equations (17) can also be expressed 
by those with respect to the new radial variables: 



d 


_ 2^ 
P 


d 
dp 


4^ 


\ 9 
tan«7r — h 
da 


cos P 
COS a 


d 

da 


sin (3 
sin a 


d ' 
dp_ 


1 


d 

W2 


_ 2^ 
P 


d 
dp 




i d 

tan a 7- 

da 


cos P 
cosa 


d 

da 


sin P 
sin a 


d ' 
dp 


1 



d f sin /9 ^ ^ ^ 

9^3 P^ [cosa da sin a dp 
Applying them to Zn,^{a,P), we obtain 

d 

tana—Zn,m{oi,P) = {2n+ \m\)Zn,m{cy, P) 

n-l 



+ 2{-ir+'{{2r+\m\ + l){2n+\m\ + l)y/^Zrm{a,P), 

r=0 

and when m > 0, 

[ cos P d sinP d \ 



cos a da sin a dp 

n— 1 



'n,±m{a,P) 

^ (-1)'^+^ {(2r + m + 2)(2n + m + 1)}^/' ^.,±(^+i)(a, P) 

r=0 

n 

+ E (-1)"^'' {(2^ + m)(2n + m + 1)}^/' /?), 



r=0 



cos/? d sin/3 d \ 

a '■ ^ r ^n>o(Qi, P) 

cos a da sma dp \ 



Y: (-ir+'-{(2r + 2)(2n + l)}^/^{z,,i(a,/3) + /3)} , 



r=0 



J sin/3 d cos/3 (9 1 
I COS a aa sm a dp \ 



(30) 



(31) 
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n— 1 



^ |(2r + m + 2)(2n + m + 1)}^/' ^.,±(^+i)(a, /?) 



r=0 

n 



^ {(2r + m)(2n + m + 1)}^/' (3), 

r=0 



. Jsin/3 a cos/5 d \ 
I cosa aa sma ap I 



n-l 



E (-ir+'-{(2r + 2)(2n + l)}^/^{z,,i(a,/3) - /?)} . (32) 

Now, expanding the radial functions ip^^ as 

oo oo 

^l\p,a,P) = Y: E <n,^(p)^n,n^(«,/3), (33) 



n=0 m=—oo 



and substituting tliem into tiie radial equations (17), we obtain a set of coupled ordinary 
differential equations for the functions Rq'^n,mip)' 

{'''^ + + + ^ " T " ^["^""^ + + ^^^^"^ + ""^^ } 

oo 

-8(£ + A) E (-l)"""'{(2n + m + l)(2r + m + l)}^/^i?^f,,±„(p) 

r=n+l 

{oo 
E (-1)"+^ {(2n + m + l)(2r + m)}'/' R2r,±^m-i)(p) 
r=n+l 

oo 1 

+ E (-1)"^^ {(2^ + m + l)(2r + m + 2)}^^ it:Jf,,±(^+i)(p) 

r=n J 

{oo 
E (-1)"+'' {(2n + m + l)(2r + m)}^/^ i?S_i),.,±(^_i)(p) 
r=n+l 

oo 

- E (-1)"^'' {(2^ + m + l)(2r + m + 2)}^/^ i?[J_i),.,±(^+i)(p) 

r=n J 
I oo 

Ti4(£ - g) E (-1)"^'' {(2^^ + m + l)(2r + m)}'/' R[i^,),Mm-i)ip) 



r=n+l 



oo 



- E (-l)"""^ {(2^ + m + l)(2r + m + 2)}^^ R\U),rMm+i)ip) 



r=n 



= 5: C(|m'Tm|)D(n,±m,r,m')pi?;i.,^,(p), 

+ P (5 + 2£ + 2A) ^ - ^ - 4[A„o + 2n(£ + A)]| i?Ji,,o(p) 
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+ A) E (-ir^{(2n+l)(2r + l)}^/^i?J;^,,o(p) 

r=n+l 

oo 

+4{i - 2g + A) E {(2n + l)(2r + 2)}'/' + 

r=n 

oo 

+z4(g - A) E (-1)"^'' {(2n + l)(2r + 2)}^/^ " 

r=n 

oo 

- g) E (-1)"+'- {(2n + l)(2r + 2)}^/^ " ^g+i),.,(-i)(p)} 

r=n 

= Ptit- E E ^^(|m'|)D(n,0,r,m')pi?^i.,^,(p), (34) 

where m > 0. 

Because the differential operators in the set of coupled ordinary differential equations 
(34) are same, it can be reduced to a system of linear algebraic equations by the spectral 
decomposition of R^^n,m{p) terms of the associated Laguerre polynomials: 

oo 

Kn,M - e-^" E &,n,n.Lf^'''-'\pl (35) 
p=0 



where 



n\ dp^ ^ ' r\\n — r)\\n — r ^ m)\ 

satisfying 



P^2^''n\p) + (m + 1 - + nL^n\p) = 0, 

= -in + 1)4+Up) + (2n + m + l)^-^^) - (n + m)L^:'\{pl 

r dpe-^p-Lt\p)L^:'\p) = (37) 

Thus, 

= -e-''^ (p+t + X + 5/2)pLf (p) . (38) 

In summary, the solution of the Schrodinger equation for the three-body problem with 
the given angular momentum £ and parity (— 1)^"*"^, where A = or 1, can be expressed 
as an infinite series: 

*a(x,y)^e-''/2E EE E /pln,„.^r^'+'Hp)^n.^(«,/5)gf(x,y), (39) 

q=\ p=0 n=0 m=—Nz 

where Nj are infinity, and the coefficients fp^q^n,m satisfy a system of linear algebraic 
equations: 

(tS^) + ^ + ^ + mPfl\,,,n,±m - [^Km + 4(^ + A)(2n + Tu) 
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+2{p + i + X + h/2f\ftln,±m + (P + ^ + A + 7/2)(p + 2^ + 2A + h)fll,,^^,^^^ 



ex 



N2 



-8(£ + A) Y: (-ir''[(2n + m + l)(2r + m + l)]^/Vp'L±. 



r=n+l 



+4(^ - 2? + A) 



N2 



Y: (-!)"+'• {(2n + m + l)(2r + m)}^^ 



r=n+l 



Af2 



+ E (-1)'^^'' {(2n + m + l)(2r + m + 2)}^/^ 



Ti^q - A) 



iV2 



E (-!)"+'• {(2n + m + l)(2r + m)}^/^ 



r=n+l 



iV2 



E (-1)"-^^ {(2n + m + l)(2r + m + 2)}'/' f^%.,),rMm+i) 



T^4(£-g) 



Af2 



r=n+l 



N2 



- E (-1)"'"' {(2^ + m + l){2r + m + 2)}'^' 
= E E C7(l"^' T m\)D{n, ±m, r, m') 

r=0 m'=-N3 

+{2p + 2£ + 2X + 5)4^i,,,^, _ (p + 2^ + 2A + S)^^^,),^,,,^,} , 

(tS^) {(^ + ^ + ^ + mPfi\,,n,o - [4Ano + 8n(^ + A) 
+2{p + £ + X + 5/2) Vp'in,o + (P + ^ + A + 7/2)(p + 2£ + 2A + 5)/,^^i,,,„,o 

Af2 

-8(£ + A) (-ir'-[(2n+l)(2r + l)]^/Vit,o 



r=n+l 



Af2 



+4(£ - 2g + A) E {(2n + l)(2r + 2)}^/^ (4^^,,,, + 

+i4(5 - A) E (-1)"-''' m + l)(2r + 2)}V^ - 

r=n 
N2 

+i4ii -q)j: {(2n + l)(2r + 2)}^/^ ' 

iV2 Afs 

^ C^(|m'|)D(n,0,r,|m'|){-p/J,^) 



r=0 m'=-N3 

+ {2p + 2i + 2X + 5)4i,,,^, - (p + 2£ + 2A + 5)/gV,),,,,,^,} , 
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where m > 0. The coefficients C{n) and D(n,m,n' ,m') were given in Eqs. (28) and 
(29). The square of the factor h'^k/ (^/2Me^) is the energy —E in the unit Me^/h^. This 
factor can be calculated by the condition that the coefficient determinant of the coupled 
linear algebraic equations vanishes. This problem of calculation can be changed to an 
eigenvalue problem. In the real calculation, the infinite series is truncated at integers Nj 
large enough to make the coefficients at the truncated terms smaller than the permitted 
error. 



V. THE HELIUM ATOM 

The problem of a helium atom is a typical three-body problem with two identical 
particles. For this problem, we enumerate the helium nucleus to be the first particle, 
and two identical electrons to be the second and the third particles. Denote by Mne and 
Me the masses of the nucleus and the electron, whose electric charges are 2e and — e, 
respectively. Thus, 

mi = 1 — 2me, m2 = = nie, 

Z\ = 2, Z2 = Zy, — —1, 

x= -<^(l-2me)/(2me) ri, y = (me/2)V2 - rg) , 

/3i = TT, /32 = -/33, (41) 

Therefore, the coefficient C{n) in Eq. (29) becomes real. Hereafter we only remain the 
leading terms and the next leading terms with respect to the mass ratio — M^/M, 
where M — Mffe + 2Me is the total mass of the hehum atom. 

C{n) J 7l/2-4(-l)"/2(i_^^/2) when n is even , , 

~ \ I (42) 

^/m~e [ -^l/2-4n(-l)('*-^)/2^e when n is odd, 

Although the spin-dependent interactions are neglected in the Schrodinger equation 
(1) for the helium atom, the symmetric effects of the spins of two electrons have to 
be taken account. The helium is called parahelium if the spinor part of wavefunction 
is antisymmetric {S — 0), and it is called orthohelium if the spinor wavefunction is 
symmetric {S — 1). Therefore, the spatial wavefunction of parahelium is symmetric with 
respect to the permutation of two electrons, and that of orthohelium is antisymmetric. 

On the other hand, in our notation, in the permutation of two electrons, y changes 
its sign and x remains invariant, namely, f3 changes its sign, p and a remain invariant, so 
that L^^"-\p) remain invariant, Z^mioi, P) changes to its complex conjugate, and Q^^'^{x,y) 
changes a factor (—1)^"^+"^. This symmetry gives some relations between the coefficients 
/p,g,n,m ^Iq. (39) such that the radial function '4>^^{p, a, (3) becomes real. 

A. S wave in parahelium 

For S wave in parahelium, the wavefunction only depends upon the radial variables, 
and is symmetric with respect to the permutation of two electrons, namely, /p^o,n,-m ~ 
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fpfi,n,m- For simplicity, we briefly denote the coefiicient /°?,,n,m by fp,n,m- 

Ni N2 N3 

*00(X, y) = e-^/^ E E E fp,n,mL^p\p) (2 - Smo) ^ [Zn,m{^, P)] ■ (43) 
p=0 n=0 m=0 

where 3fJ[Z] denotes the real part of Z, and S^q is the Kronecker 5 function. From Eq. 
(40), the coefficients fp^n,m satisfy a real system of linear algebraic equations with the 
equation number {Ni + 1){N2 + 1)(A^3 + 1): 



y/2rn'eMe'^ 



,n,m ( 



2" 

fp,'' 



N2 N3 r, r 



E E [C(|m^ — m|).D(n, m, r, m^) + C(m^ + m).D(n, m, r, — m^)] 

r=0 m'=0 



X 



{-P/(P-I),r,m' + (2p + h)fp^r,m' " (P + 5)/(p+i),^,^' } . (44) 

Therefore, the wavefunction \E'oo(x, y) is real. 

The square of the factor Ti^k/ {^/2mlMe^) in Eq. (44) is the binding energy —E in the 
unit e^/ao, where ao = ti^/{M^e'^) is the Bohr radius. This factor can be calculated from 
the condition that the coefficient determinant of Eq. (44) vanishes. The coefficient matrix 
on the left hand side of Eq. (44) is a direct product of a (A^i + l)-dimcnsional matrix and 
a unit matrix. Removing it to the right hand side by right-multiplying with its inverse 
matrix, we reduce the problem of calculating energy into an eigenvalue problem. In the 
numerical calculation, we take A'^i = 7, A^2 = 8 and A^3 = 16, and obtain the energy for 
the ground state of a helium atom: 

-E^^''" = 2.88935 (eVoo) = 78.62eV, (45) 

where the physical constants e, h, c, Mg, and are quoted from the particle physics 
booklet, 1998. The experimental value [2] for the minimum energy required to remove 
both electrons from a helium atom is 2.90351e^/ao ~ 79.00eV. The relative deviation of 
the calculated energy from the observed value is in the thousandths. 

We also obtain the eigenfunction for the ground state of helium, expressed as the series 
(43), which is normahzed such that the largest coefficient in modulus is one (/o,o,o)- Those 
coefficients whose absolute values are larger than 0.001 are listed as follows: 

/o,o,o = 1.0000, /o,o,i = 0.0411, /o,o,2 = -0.1047, /o,o,4 = 0.0209, /o,o,6 = -0.0076, 

/o,o,8 = 0.0030, /o,o,io = -0.0016, /o,i,o = -0.0899, /o,i,i = -0.0032, /o,i,2 = 0.0227, 
/o,i,4 = -0.0069, /o,i,6 = 0.0032, /o,i,8 = -0.0015, /o,2,o = 0.0204, /o,2,2 = -0.0075, 
/o,2,4 = 0.0029, /o,2,6 = -0.0016, /o,3,o = -0.0069, /o,3,2 = 0.0032, /o,3,4 = -0.0015, 
/o,4,o = 0.0029, /o,4,2 = -0.0016, /o,5,o = -0.0015, /i,o,o = -0.0440, /i,o,i = -0.0058, 
/i,o,2 = 0.0293, /i,o,4 = -0.0082, /i,o,6 = 0.0031, /i,o,8 = -0.0014, /i,i,o = 0.0271, 
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/i,i,2 = -0.0085, /i,i,4 = 0.0030, /i,i,6 = -0.0014, /i,2,o = -0.0081, /i,2,2 = 0.0031, 

/i,2,4 = -0.0014, /i,3,o = 0.0030, /i,3,2 = -0.0014, /i,4,o = -0.0013, /2,o,o = 0.0036, 

/2A2 = -0.0037, /2,o,4 = 0.0020, /2,i,o - -0.0039, /2,i,2 = 0.0019, /2,2,o = 0.0019. 

(46) 

The remaining coefficients can be obtained from us upon request. 

Changing the electric charge of the nucleus, we can obtain the energies of the ground 
states of the helium-like ions. The calculation results and the observed results are listed 
in Table 1. 



Table 1 Calculated and observed energies of ground states 
for helium-like ions (in the unit e^/ao) 



Ion 


Calculated " 


Observed ^ 


Relative error 


H- 


0.520618 






He 


2.88935 


2.90351±0. 00004 


0.005 


Li+ 


7.25064 


7.27980±0. 00050 


0.004 


Be++ 


13.6055 


13.65600±0.00100 


0.004 


B+++ 


21.9543 


22.03200±0.00150 


0.004 


C++++ 


32.2973 


32.40700±0.00400 


0.003 



a. By Eq. (43) with A^i = 7, A^2 = 8 and = 17. 

b. See Refs. [22] and [23] 

C. {Eobs — Ecai)/Eobs 



B. S wave in orthohelium 



For S wave in orthohelium, the wavefunction also depends upon the radial variables, 
but is antisymmetric with respect to the permutation of two electrons, namely, /° q „ -m — 
-fpl,n,m- Foi^ simphcity, we briefly denote the coefficient /°?,,n,m by -igp,n,m- 



*oo(x, y) = 2e-^/' E E E 9p,n,mL^\p)^ [Zn,m{^, /?)] , 

p=0 n=0 m=l 



(47) 



where '<S[Z] denotes the imaginary part of Z. From Eq. (40), the coefficients gp^n,m satisfy 
a real system of linear algebraic equations with the equation number {N^ + 1)(A^2 + l)-^3: 



P + 2 ) P9p-l,n,m 



9p,n,r: 



,n,m ( 



N2 N3 

E E [C{\m' — m\)D{n,m,r,m') — C{m' + m)D{n,m,r, —m')] / ^yn^e 

r=0 m'=l 

X { + (2p + ^)9p,r,m' - (P + 5)5((p+i),^,^/} . 



(48) 
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Therefore, the wavefunction ^oo(x, y) is also reaL Similarly, the energy of the ground 
state of the orthohelium can be calculated as an eigenvalue problem. In the numerical 
calculation, we take Ni = 7 , N2 — 8 and = 17, and obtain the energy for the ground 
state of the orthohelium: 

-E^^''" = 2.08039 (eVao) = 56.61eV. (49) 

The experimental value for it is —Eqq^ = 2.17524 (e^/ao)=59.19eV. The relative error is 
4%. We will discuss this problem in Sec. VI. 



C. P wave in parahelium 



There are two sets of solutions for the P wave with the different parities. The eigen- 
functions Q^^ of the angular momentum [see Eq. (14)] are written explicitly as follows, 
where g = or 1 when A = 0, and q = I when A = 1: 



y)^xi + ix2, Qo°(x, y) = yi + m, 

QVi^, y) = (xi + 1x2) ys - xs {yi + iy2) ■ 



(50) 



From the condition that the spatial wavefunction in parahelium is symmetric with respect 
to the permutation of two electrons, we have 

f = fio _ flO 
Jp,n,m — Jp,l,n,m J p,l,n,—mi 

'■yp,n,m — Jpfi,n,m ~ J p,0,n,-m-> 
~ihp^n,m = fp,l,n,m ~ ~fp,l,n,-m- (^1) 

Thus, the wavefunctions are 

Ni N2 N3 

*io(x, y) = E E E fp,n,mL^p^ ip) (2 - S^o) ^ [Zn,m{c^, P)] (x, y) 

p=0 n=0 m=0 

+ 2e-''/2 E E E 9p,n,mLf^ (P)» [^n,n^(«, P)] QI'{^, y), 
p=0 n=0 m=l 

N7 Na JV.) 

*n(x,y) =2e-''/2EE E Km(«, /?)] g}^(x, y). (52) 

p=0 n=0 m=l 

Prom Eq. (40), the coefficients satisfy a real system of linear algebraic equations, so that 
the coefficients are real: 



,^/2^Me2){(^+2j^^^-''"'^ 



4A„^ + 4(2n + m)+2(p+^ 
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N2 



fp,n,r, 



+ (p + 7)W,m-8 E (-l)"+M(2n + m + l)(2r + m + l)]^/Vp,r.,m 

r=n+l 



16 



-4 



N2 



(1-M E (-ir+''{(2n + m + l)(2r + m)}^/Vp,r,m-i 



r=n+l 



N2 

+ (1 + <5mo) E (-1)"^' {(2^ + m + l)(2r + m + 2)}^/' 



-4 



E (-ir+'{(2n + m+l)(2r + m)}'/'^p,,,^_i 

r=n+l 



- (1 + (5„o) E (-l)"^' {(2^ + m + l)(2r + m + 2)}^/' 



N2 N3 9 _ r 



r=0 m'=0 



{C(|m' — m\)D{n, m, r, m!) + C{m' + m)D{n, m, r, —m')} 



X {-p/p — l,r,m' + (2p + 7)/p,r,„/ - (p + 7)/p+i,r,„j/} , 



4A„^ + 4(2n + m) +2 (p 



9p,r 
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+ P + (P + 7)gp+i,n,m - 8 E (-1)"^' [(2^ + m + l)(2r + m + 1)]^/' ^p,,,^ 



+ 4 



E (-1)"+^ {(2n + m + l)(2r + m)}^/' 



r=n+l 



+ E (-1)"^' {(2n + m + l)(2r + m + 2)}^/' ^^,,,^+1 



+ 4 



iV2 



N2 



E (-l)'*^^ {(2n + m + l)(2r + m)}^/' 



r=n+l 



E (-1)"+'' {(2n + m + l)(2r + m + 2)}^/' /p,,,^+i 



Nr, Ne 

E E — ?7T.|)D(n, m, r, m') — (7(m' + m)D(n, m, r, — m')} /y^TTT^ 

r=0 m'=l 

X {-P5'p-l,r,m' + (2p + 7)gp^r,m' " (p + 7)5f(p+i),r,m' } , (53) 
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^V2^Me^)U^+2)^^-^- 



/ 9 

4A„^ + 8(2n + m) + 2 ( p + - 



/ 11\ 1 

+ P+V (^' + 9)Vi,n,m-16 E (-ir+M(2n + m + l)(2r + m + l)]^/' V,™^ 

^ ^ ^ r-=n+l J 



TVs iV') 



= E E {<^(|^' — ''7T'|)-D(n, m, r, m') — (7(m' + m)D(n, m, r, — m')} /y^TTT^ 



=0 m'=l 
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X {-p/lp_l,r,m' + (2p + Q)hp^r,m' " (P + ^)hp+l,r,m'} , (54) 

The energy can be calculated as an eigenvalue problem. In the numerical calculation, we 
take Ni ^ N4 = Nt = 7, N2 = = Ns = 8, iVg = Ne ^ Ng = 17, and obtain the 
lowest energies of the P wave with the odd parity and the even parity for the parahelium, 
respectively: 

-£;fo''" = 2.02095 (eVao) = 54.99eV, (58) 

= 0.581291 (eVao) = 15.82eV. (59) 

In the shell model of the atomic physics, the electrons in a helium atom are supposed 
to fill in the energy levels of a hydrogen-like atom with the electric charge +2e according 
to the exclusion principle. In this model, the lowest-energy state of the P wave with the 
odd parity in parahelium is explained as the compound state of one 15" electron and one 
2P electron. The observed energy is -E^^^ = 2.1238(e^/ao) =57.79eV. The relative error 
for the calculated energy —E^q^'^ is 5 %. 

On the other hand, it is worthy to pay more attention to the existence of the lowest- 
energy state of the P wave with the even parity in the parahelium. This state is forbidden 
for the electric dipole transition to the ground state of the parahelium, but it is allowed to 
the P wave state with the odd parity in the parahelium. The energy difference between 
two P wave states with different parities is 

A£; = 1.43966 (eVao) = 39.18eV = 315969cm-\ (60) 

It is larger than the ionization energy (24.58eV) of the helium atom, which is the minimum 
energy required to remove one electron from the ground state of a helium atom. In the shell 
model, this P-wave state has to be composed of two P-wave electrons. According to the 
Clebsch-Gordan coefficients [20] , the angular part of its wavefunction is antisymmetric [see 
Eq. (50)]. Therefore, this state should be explained in the shell model as the compound 
state of one 2P and one 3P electrons, because its radial function has to be antisymmetric. 
However, the binding energy —E of the compound state seems to be smaller than our 
calculated value. We will further discuss this state in Sec. VI. 

D. P wave in orthohelium 

We sketch the calculation on the lowest-energy states of the P wave in orthohelium. 
Prom the condition that the spatial wavefunction in the orthohelium is antisymmetric 
with respect to the permutation of two electrons, we have 



^fp,n,m 


= fio 

— J p,l,n,m 


_ rlO 

Jp,l,n,— 


9p,n,m 


= fio 

— Jp,0,n,m 


_ rlO 

J p,0,n,— m' 




= f^^ 

J p,l,n,m 


- f^^ 

J p,l,n,—m' 



(61) 



Thus, the wavefunctions are 

^rio(x, y) = 26-"/' E E E fp,n,mLf\p)Q [Z^,m{a, /?)] gi°(x, y) 

p=0 n=0 m=l 
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+ e-'/' E E E 9p,n,mL^p\p) (2 - S^o) ^ [Zn,m{<^, Z^)] Qo°(x, y), 

p=0 n=0 m=0 

N7 Ns Ng 

*n(x, y) = E E E K,n,mL^p\p) (2 - <^mo) 3? [^n,m(«, /?)] gi^(x, y). (62) 

p=0 n=0 m=0 

From Eq. (40), the coefficients satisfy the real systems of hnear algebraic equations, so 
that the coefficients are real. The real systems of hnear algebraic equations are similar to 
that for the parahelium, except for changing some signs. We will not list those equations 
here, but list the calculated results for the energy, where the series are truncated at 
Ni = N4 = NT = 7, N2 = N5 = Ns = 8, = Nq = Ng = 17, 

_Eortho ^ 2.04388 (eVao) = 55.62eV, (63) 

_^ortho ^ 0.710413 (eVao) = 19.33eV. (64) 

In the shell model of the atomic physics, the lowest-energy state of the P wave with 
the odd parity in orthohelium is also explained as the compound state of one 15* electron 
and one 2P electron, and that with the even parity is explained as the compound state of 
two 2P electrons. The observed energy is -E^^' = 2.1332(eVao)=58.05eV. The relative 
error for the calculated energy —E'[q^"' is 4 %. 



VI. DISCUSSIONS 



In this paper we have presented an improved hyperspherical harmonic method for 
directly solving the quantum three-body problem. Applying this method to the problem 
of the helium atom, we have calculated the lowest energies and the wavef unctions of some 
states of the S wave and the P wave in the parahelium and in the orthohelium. The 
relative error for the energy of the ground state of the parahelium is in the thousandths, 
but the remaining results have the relative error of a few percents. 

The main approximation, which we have made in our calculation, is to take the finite 
Nj in the series for the wavef unctions. If the coefficients of the truncated terms in the series 
are much smaller than the needed precision, the calculated results will be satisfactory. For 
example, for the ground state of the helium, when Ni = 7, N2 = 8 and — 16, we find 
that all the coefficients at the truncated terms are much less than 0.001: 



/7,o,o = 1.046 X 10-^ /o,8,o = 2.995 x 10-^ /o,o,i6 = 2.983 x 10" 



(65) 



Therefore, the relative error of the calculated energy must be in the thousandths. It also 
can be seen in the calculated results of the energy when different Nj are taken: 



Note that, when 



2.88331(eVao), 
2.88805(eVao), 

2.88935(eVao), 
2.88935(eVao), 
2.88972(eVao), 

= 8 and Nj, = 16, the calculated energy when Ni 



when 


Ni 


= 5, 


N2 


= 4, 


A^3 


= 8, 


when 


Ni 


= 5, 


N2 


= 6, 




= 12, 


when 


Ni 


= 5, 


N2 


= 8, 


N3 


= 16, 


when 


iVi 


= 7, 


N2 


= 8, 


Ns 


= 16, 


when 


Ni 


= 7, 


N2 


= 8, 


Ns 


= 25. 



(66) 



5 is the same as 



that when Ni = 7. From those data we may conclude that the finite Nj we have taken 
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for the calculation of the energy of the ground state of the helium atom are suitable for 
the relative precision in the thousandths. 

From the viewpoint of calculation, to raise A^^i and will only increases the calculating 
quantity, but when N2 raises, the formula for calculating the coefficients D{n,m,n' ,m') 
becomes much complicated. Now, it is hard for us to raise N2 more than eight. It is the 
reason that we have to accept the precision of a few percents for the remaining calculated 
energies, and to leave the improvement of the precision in the future. For example, for the 
lowest-energy state of the S wave in the orthohelium, when Ni = 7 , N2 = S and = 17, 
the coefficients at the truncated terms are 



^7,0,1 = 3.430 X 10-^ ^0,8,1 = 1-663 x 10-^ ^o,o,i7 = 1-650 x 10" 
and the calculated energies are as follows when different Nj are taken: 



(67) 



lortho 
00 



— < 



' 2.06704(eVao), 
2.06762(eVao), 
< 2.07932(eVao), 
2.08039(eVao), 
2.08396(eVao), 


when 


iVi 


= 5, 


N2 = 6, 


N3 


= 13, 


when 


iVi 


= r. 


N2 = 6, 


Ns 


= 13, 


when 


N, 


= 5, 


N2 = 8, 


Ns 


= 17, 


when 


Ni 


= 7, 


A^2 = 8, 


Ns 


= 17, 


when 


Ni 


= 7, 


N2^8, 


N3 


= 25. 



(68) 



The relative error is a few percents. 

We would like to pay more attention to the lowest-energy state of the P wave with 
the even parity in the parahelium. In order to show the property of the coefficient /ip,n,m 
in the series (52) as subscripts increase, we list some coefficients /ip,n,m when A^i = 7, 
N2 = 8, and N3 =17: 



^0,0,1 — 1, ^1,0,1 — ~ 


0.4015, 


^2,0,1 — 


0.1362, /i3,o,i = - 


0.0446, 


^4,0,1 = 0.0153, /i5,o,i = 


-0.0055, 


^6,0,1 


= 0.0019, 


^7,0,1 = 


: -0.0006, 


^0,1,1 = -0.3583, /io,2,i 


= 0.1460, 


^0,3,1 


= -0.0692, 


^0,4,1 


= 0.0374, 


ho,5,i = -0.0227, /io,6,i 


= 0.0152, 


^0,7,1 


= -0.0117, 


^0,8,1 


= 0.0086, 


^0,0,2 = 0.0501, /io,o,3 = 


: -0.3675, 


^0,0,4 


= -0.0104, 


^0,0,5 


= 0.1460, 


^0,0,6 = 0.0038, ho,o,7 = 


= -0.0679 


^0,0,8 


= -0.0013, 


^0,0,9 


= 0.0352, 


^0,0,10 = 0.0006, /io,o,ii = 


-0.0201, 


/io,o,i2 = -0.0003 


/io,o,i3 = 0.0122 


^0,14 = 0.0002, /io,o,i5 = - 


-0.0078, 


/i0,0,16 - 


= -0.0001, 


^0,0,17 


= 0.0051, 



We also list some calculated results for the energy when N2 — 8: 



rppara 



0.576911 
0.579576 
0.580345 
0.577204 
0.580270 
0.581291 
0.581840 



e^/oo) when A^i 

e^/flo) when A^i 

e^/ao) when Ni 

e^/ao) when A^i 

e^/ao) when A^i 

e^/ao) when A^i 

e^/ao) when A^i 



5, and A^3 

5, and N3 

5, and N3 

7, and A^3 

7, and A^s 

7, and A^a 

7, and N3 



9 

13 
17 
9 

13 
17 
25. 



(69) 



(70) 
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Due to the restriction in N2, the relative error of the calculated energy is in the hundredths. 
However, the existence of a P wave state of even parity in the parahelium with the energy 
—E ~ 0.58(e^/ao) is reliable. Note that this state is forbidden for the electric dipole 
transition to the ground state, but it is allowed to the P wave state with odd parity. 
The energy difference in the allowed transition is around 39.18eV. We are waiting for the 
observation of this transition in experiments. 

It seems that the existence of this state conflicts with the shell model in the atomic 
physics. In the shell model two electrons in the helium atom arc supposed to fill in the 
energy levels of a hydrogen-like atom according to the exclusion principle. The lowest- 
energy state of the P wave with the even parity has to be composed of two P wave 
electrons. From the shell model, two 2P electrons cannot compose a P wave state with 
symmetric spatial wavefunction, and this state has to be explained as the compound state 
of one 2P and one 3P electrons. In this state each electron moves in the electric field of 
the nucleus, screened by the other electron: 

2 • 22 • ao 2 • 32 • ao ■ 

If one electron moves completely inside of the other electron, Zi = 2 and = 1, or 
vice versa. In the real case, Zi is less than 2 and Z2 is larger than 1. We assume that 
Zi = 2 — r, Z2 = 1 + r, and < r < 1. Thus, the energy —£^10 = —E{t) is a function of 
the parameter r, and given in the unit e^/ao as follows: 

-E(Q) = 0.556, - E{Q.l) = 0.518, - E{Q.2) = 0.485, - E{Q.2,) = 0.455, 

-E(0.4) = 0.429, -E(0.5) = 0.406, - £(0.6) = 0.387, - £(0.7) = 0.372, 

-£(0.8) = 0.360, -£(0.9) = 0.352, - £(1) = 0.347. 

This estimation is very rough, but it shows qualitatively that this P wave state seems not 
to be the compound state of one 2P and one 3P electrons in the shell model. If it is true, 
the shell model in the atomic physics is challenged even qualitatively. We sincerely hope 
that this P wave state with the even parity in the parahelium will be studied further in 
both theoretical and experimental physics. 
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Appendix PROOF OF SOME IDENTIFIES 

In this appendix we will prove some identities for the alternating series used in this 
paper. 

V (-ir(^ + ^ + ^)! ^ ._l^n 

h r\{n-r)\{m + r)\ ^ ' ' ^ ' 

Proof: From Newton's binomial expansion 

oo 

(l + a;)-"= E 



Q m!(n — m)!' 

-l)™(n + m- l)!,x'" 
m=o m!(n-l)! 
we have 

oo 

(-i)V = (1 + x)-^ = (1 + x)"-"-^ 



t=0 



r=o ri[n-ry. s\n\ 



22 



to^ ' h r\{n-T)\{t-r)V 

where the summation index s is replaced hy t — r-\-s. Comparing two sides of the equahty, 
we obtain 

y (-ir(n + t-r)! 
^ r\{n — r)\{t — r)\ 

Letting t = n + m, and then replacing the summation index r by n — r, we proved Eq. 
(Al). Q.E.D. 

^(„, ,) , f (-l);(« + m + r)! _ (-!)■■(* - »)- , (^2) 

where (a)„ = a(a + 1) ■ ■ ■ (a + n — 1) and (a)o = 1. 

Proof: Prove Eq. (A2) by induction. It is obvious that Eq. (A2) holds when n = 
and n — 1. Now, assuming that Eq. (A2) holds for n < s, we will prove it holds when 
n — s-\-l. In fact, 



F(.+ l,™,i)(™ + ( + .+ l) = i; (-l)%' + l+™ + '-)!(m + t + r) + (.+ l-r) 

r= 



r!(s + 1 — r)!(m + r)! (m + t + r) 



= (_l)^+i + V (-l)'(g + m + ry.{m + t + r) + {s + l-t) 
^ r\{s — r)\{m + r)\ (m + t + r) 

= -F(s,m,i)(i-s- 1). 

Removing the factor (m + t + s + l) from the left to the right, we proved that Eq. (A2) 
holds for n = s + 1. Q.E.D. 



^ r!(n — r)!(m + r)! {m + t + r)^+i 



where 



[m + t)n+£+l 

p(0) _ -I 

=t(m + t), 
Pi^^^ = -t + n{n + m + l), 
pP ^{t + l)t{m + t){m + t + l), 
Pf ^ = 2t{m + t)[-{t + l)+n{n + m + 1)], 
^2^2) = + 1 _ 2n(n + m + 1)] + (n - l)n(n + m + 2)(n + m + l)/2, 
P^^^ ^{t + 2){t + l)t{m + t){m + t + l)(m + t + 2), 
pf ) = 3(i + l)i(m + i)(m + i + l)[-{t + 2) + n(n + m + 1)], 
P^^^ = 3t(m + + + 2 - 2n(n + m + 1)] 
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+ n{n - l)(n + m + l)(n + m + 2)/2}, 

pP = t{t + l)[-{t + 2) + 3n{n + m + l)] - 3tn{n - l){n + m + l){n + m + 2)/2 

+ n{n — l)(n — 2)(n + m + l)(n + m + 2)(n + m + 3)/6, 

Po^^^ = {t + 3){t + 2){t + l)t{m + t){m + t + l){m + t + 2){m + t + 3), 

) = 4(t + 2)(t + l)t(m + t)(m + i + l)(m + t + 2)[-{t + 3)+n{n + m + 1)], 

P2(^) = 6{t + l)t{m + t){m + t + l){{t + 2)[t + 3- 2n{n + m + 1)] 

+ n{n- l)(n + m + l)(n + m + 2)/2}, 

P^^^ = 4t(m + + + 2)[-(t + 3) + 3n{n + m + 1)] 

-3(i + l)n(n- l)(n + m + l)(n + m + 2)/2 

+ n(n - l)(n - 2)(n + m + l)(n + m + 2)(n + m + 3)/6}, 

p]^) = f(i + l){t + 2)[i + 3 - 4n(n + m + 1)] 

+ 3t{t + l)n(n - l)(n + m + l)(n + m + 2) 

- 2tn{n - l){n - 2){n + m + l)(n + m + 2)(n + m + 3)/3 

+ n(n - l)(n - 2)(n - 3)(n + m + l)(n + m + 2)(n + m + 3)(n + m + 4)/24, 

The remaining coefficients can be obtained from us upon request. 

Proof: When £ = 0, G(n, m, t, 6, 0) = F(n, m, t). When £ = 1, we have 

G(n, m, t, t, 1) = (-ir(^-^ + l)n-i p(i) = ^(^^ ^ + 

(^m + tJn+2 

m, M + 1, 1) = ^"'7:"/'^"-^ {P^'' - Pi'' } = F{n, m, t). 
(m + 1)„+2 

Pi^^^ and P^^^ can be solved from those relations. When i = 2, we have 

G(n, m, t, t, 2) = (-in^- ^ + 2)n-2 (2) ^ ^^^^ ^ ^ 3), 

(m + i)„+3 

G(n, m, t, t - 1, 2) = ^'^^''^^ ' " ^ "^"'^ [P^^^ + Pf } = G(n, m, t + 1, t - 1, 1), 
G(n, m, M + 1, 2) = (-^^(^ " ^ + 2)n-2 |p(2) _ ^(2) ^ ^pC^) j ^ ^ 
Gin, m,t,t + 2, 2) = i-^t -n-^^)n-. r (2) _ ^^(2) ^ gp(2)} 

= G{n, rrijtjt + 3, 1). 

The first three relations are used for calculating the coefficients P2^^ , P^^^ and Pq^^ , and the 
last one is for check. Similarly, the coefficients Pp can be calculated from G{n, m, t, h, i'), 
where £' < £. 
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